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1. Overview 



An important contribution of lattice QCD to phenomenology is the calculation of quark masses. 
Here, we discuss the mass calculation of the heavy quarks bottom and charm. Our method com- 
bines Monte Carlo calculations of heavy-light mesons with lattice perturbation theory. This first 
section provides an overview of the calculation which uses the Fermilab method for heavy quarks. 
Sections 2 and 3 review the non-perturbative and perturbative aspects, respectively. We conclude 
with preliminaiy results for the bottom quark mass. 

Because the lattice violates Lorentz (Euclidean) invariance, the energy-momentum relation- 
ship, 

E\p)=ml + '^p^ + ..., (1.1) 

is distorted [1]. The mass mi is called the rest mass and m2 the kinetic mass. They are defined, for 
quarks or hadrons, as [1] 

mi=E{0) „^2= . (1.2) 



p=0 



Both m\ and m2 provide a means for determining the quark mass. The two methods must yield the 
same result in the continuum limit, providing a cross-check on our results. 

The first method uses the meson binding energy to arrive at the quark mass. For the Fermilab 
method, it has been shown that [2], 

Ml - mi = Mexpt - mpoie mpoie = nil + {M^xpt -Mi) (1.3) 

where Mi is the heavy-light meson rest mass calculated on the lattice, Mexpt is the experimentally 
measured meson mass, mi is the lattice heavy-quai^k pole rest mass and mpoie is the continuum 
quark pole mass. This equation holds up to discretization enors in operators of dimension six 
and higher, and to truncation eiTor in mi, when that is defined perturbatively. Specifically, the 
leading mismatch in the Lagrangians comes from the hyperfine interaction, the Darwin term and 
the spin-orbit interaction [2]. Equation (1.3) also holds for spin-averaged mesons. Mi and Mgxpt, 
so spin-dependent discretization effects can be eliminated by using them. This leaves the Darwin 
term as the leading source of discretization errors. Calculating the left-hand side of Eq. (1.3) on the 
lattice and using PDG [3] values for Mexpt, we arrive at a value for mpoie. We refer to this method 
as the rest-mass method. 

Alternatively, we can use the kinetic mass which the Fermilab method identifies with the 
physical quark mass, am2 = ampoie. We use the ratio Mexpt/«M2 to set the (inverse) lattice spacing. 
The quark mass is then 

Mexpt am2 

mpoie = am = -=- Mexpt- (1-4) 

aM2 aM2 

The advantage of setting the lattice spacing this way can be seen in the second equality. It shows 
how a mistuning of the heavy quark mass cancels in the ratio am2/aM2. We refer to Eq. (1.4) as 
the kinetic-mass method. 
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Figure 1: Spin-averaged meson rest mass versus the light valence mass in units of rj [8] on the 0.09 fm 
lattice. The straight line on top of the three heaviest data points demonstrates clearly the linearity of the data 
near the physical strange quark mass. 

2. Non-perturbative Elements 

We use the MILC 2+1 flavor lattices [4], which have asqtad sea quarks [5] and improved 
gluons [6]. The calculation includes three lattice spacings of approximately 0.09, 0.12, and 0.15 
fm. The ratios of the nominal up-down to strange quark masses range from niuj/ms = 0.1 to 0.3 
or 0.4 depending on the lattice spacing. The heavy bottom and charm quarks are simulated using 
the Fermilab method [1]. An asqtad quark is used for the meson's light valence quark. 

Meson rest masses are calculated using constrained curve fits to two-point correlators [7]. We 
determine masses of both the pseudoscalar and vector mesons (e.g., and B*) and spin average. 
To determine the kinetic meson mass, M2, we first determine the spin-averaged value of the energy 
at several values of momenta and then fit to the dispersion relation, Eq. (1.1) for mesons. 

Heavy-light mesons used in this calculation have a strange valence quark. This allows us 
to avoid a chiral extrapolation in the valence mass. Still, simulations are not done exactly at the 
physical strange mass. To reach it, we linearly interpolate between two neighboring points, which 
is validated by Fig. 1 . 

3. Perturbation Theory 

We use one-loop perturbation theory results to obtain the quark pole masses mi and m2 from 
the bare mass [9, 10]. We use the V-scheme [11, 12] for the strong coupling a{q*) = g^{q*) /An 
and determine its value in the manner described by Mason et al. [13]. The scale q* should be 
chosen to be the typical momentum of a gluon in the loop. To determine that momentum, we use 
the method introduced by Brodsky, Lepage and Mackenzie (BLM) [11, 12] extended to include 
cases where the one-loop contribution is anomalously small [14]. BLM defines q* by 



\nq 



,*2 



jdUf{q)Hq^) 
Jd'qfiq) 



(3.1) 
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where f{q) is the one-loop integrand. In cases where the one-loop integral is anomalously small, 
or zero, this expression is inappropriate. If the one-loop integral is zero. Reference [14] defines q* 

. (3 2) 

Solutions for small but non-zero one-loop integrals allow for a continuous transition from Eq. (3.2) 
toEq. (3.1). 

It is known that the pole mass is not a good choice of scheme due to renormalon ambigui- 
ties [15]. Therefore, we use a short-distance mass which is designed to run sensibly at low renor- 
malization scales [16]. We use the potential subtracted mass [17] because it is additive and therefore 
works well with the rest-mass method as shown below. It is based on the static quark potential and 
introduces a separation scale /ij. At the one-loop order 

mps(At/) = mpoie - + (3.3) 

71 47r 

where Cf = 4/3 and Aqcd < M/ < "Vark- 

Generically, the expansion for the quark rest mass from the lattice can be written as 

[0] I 2 [1] , 4 [2] , .X 

am\ = am\ -\-g am\ -\-g am\ + • • • (3.4) 
and similarly for ni2- For the rest mass we can define 

amiMl^f) = amf^+g'^{q*) jam^ - ^^^^^ | , (3.5) 
such that Eq. (1.3) yields 

mps(At/) =mi,ps(At/) + (Mexpt-Mi), (3.6) 
where mps(/X/) is the continuum quark mass in the potential subtracted scheme. For m2 

Mexpt Cp^fg^jq*) 
mps(M/)=«-,-=--— — . (3.7) 

Below, our final quark masses for bottom are quoted at the scale jU/ = 2 GeV. We reach this 
scale in two ways, which treat higher order effects differently. The first uses a fixed jXf and the 
calculation of q* and g^{q*), as described above. We call this the q* method. The second has /i^ 

chosen such that the one-loop coiTcction is zero, e.g. for the rest-mass method jj.'^ = An^n^l^ /Cp- 
The resulting mps(/i^) is then run to the final /ij by using the two-loop renormalization group 
equation for the PS mass: 

mps{iif)-mps{ii'f) = — [At}o;(At/){l+2jSoa(At})}-At/«(At/){l+2j8oo;(At/)}] , (3.8) 
where jSo = 4^(11 — 2/3 «/). We call this the zero-and-run method. 
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Figure 2: (a) Results for the bottom quark mass in the potential subtracted scheme for three lattice spacings. 
Offsets on the x-axis are for clarity, (b) Spin-averaged (rest) mass at two values of the valence mass and 
three different values of sea mass ratios mu.d/ms. Data are from the 0. 12 fm lattice. Offset on the .jc-axis is 
for clarity 



4. Preliminary Results 

Preliminary results for the bottom quark mass were presented last year [19]. Several improve- 
ments have been made since then. First, we now include the 0.09 fm lattice and omit the 0.18 fm 
lattice. Second, we have added the kinetic-mass method. Third, we use improved scale setting for 
(x{q*) [14], and finally we use the zero-and-run method for the mass scale. 

Figure 2(a) shows updated results for the bottom quark mass. Error bars include uncertainties 
from statistics, chiral sea-quark effects and the determination of the lattice spacing. Most noticeable 
are the large error bars on the kinetic mass results (open symbols) versus the rest-mass results (filled 
symbols). These are primarily statistical; we hope to reduce them in the future so that the kinetic- 
mass method is a stronger cross-check of the rest-mass method. 

Table 1 lists the percent uncertainties in the calculation of the bottom quark mass from the 
rest-mass method. The uncertainty due to the truncation of the QCD perturbation theory clearly 
dominates. To estimate this, we take the spread in results for the quark mass from the q* and 
zero-and-run method on the 0.09 fm lattice. This yields a 4% uncertainty which is consistent with 

We expect the dependence on the mass of up-down sea quarks to be mild for heavy-strange 
meson masses. Fig. 2 (b) shows a typical set of meson (rest) masses for three different sea-quark 
ensembles. The up-down quark mass gets smaller from the blue triangle to the green square to 
the red circle. For our central value, we use the meson mass from the ensemble with the smallest 
up-down quark mass and we take the largest spread in values as an estimate of the uncertainty due 
to sea-quark effects. For example, in the plot shown we use the red circle for our central value and 
the difference between it and the green square as an estimate of the effect of m^^d 7^ 0- 

For the heavy-quark discretization error we consider only the contribution from the Darwin 
term, since spin- averaging removes the hyperfine and spin-orbit interactions. The coefficient for 
this term from both the Sheikoleslami-Wohlert and continuum actions is known. Using the dif- 
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Source 


percent error 


statistical 


0.1 


lattice spacing determination 


0.4 


heavy-quark tuning 


0.5 


sea quark effects 


0.7 


strange mass tuning 


0.2 


perturbation theory truncation 


4 


light quarks and glue 


1 


heavy quark discretization 


0.6 


total 


4.3 



Table 1: Percent uncertainties in the bottom quark mass. These are added in quadrature to arrive at the total. 



ference, fg, one can estimate the uncertainty as a^/fAgf^p given the lattice spacing a and an 
estimate for Aqcd [18]- To get an estimate for Aqcd> we fit two subsets of the rest-mass results 
(0.09 and 0.12 fm; 0.09 and 0.15 fm) to an 0{a^) ansatz. The average of these results gives 
Aqcd = 1-3 GeV. Evaluating fE^^QQ^ at the 0.09 fm spacing then yields a 0.6 percent error. 
Although Aqcd = 1-3 GeV is high, the resulting error is small and so we conservatively take this 
as the uncertainty due to heavy-quark discretization. 

The lattice spacing (or r\ [8]) determination and strange mass tuning are done by the MILC 
Collaboration [20]. We use r\ = 0.318(7) fm. We also include errors due to an estimated 10% un- 
certainty in the tuning of and an 8% mistuning of the heavy quark. We estimate the uncertainty 
due to discretization of the light quarks and gluons as a a^I^QQYy and o'Kqqy) with Aqcd = 1 GeV 
and quote the larger of the two. 

From Fig. 2 (a), the rest-mass method clearly has much smaller errors. In addition, the lattice 
spacing dependence is see to be mild. For these reasons, we take the average of the two results from 
the q* and zero-and-run methods for the rest-mass on the 0.09 fm lattice as our central value for 
the bottom quark mass. Uncertainties in Table 1 are added in quadrature to arrive at the total. We 
have then for the bottom quark mass in the potential subtracted scheme mb,ps(2 GeV) = 4.32(19) 
GeV. For comparison, a QCD sum rule calculation [21] obtains mh,ps{2 GeV) = 4.52(6) GeV and 
mh{mh) = 4. 19(6) GeV in the MS scheme. 
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